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The dependence of the Josephson plasma frequency ui p in Bi2Sr2CaCu2 0s on a tilted field H is 
reported. Measurements over a large range of B and tilt angle 9 allow a detailed comparison with a 
recent calculation by Koshelev. With a slight modification of the model, close agreement is obtained. 
From the fits, we find values for the in-plane correlation length and the zero-field critical current 
density J c o (4,600 A/cm 2 at 30 K). An analogy to Bragg diffraction is described, as well as a picture 
for the fractional-exponent behavior of u v versus H. 



PACS numbers: 74.72.Hs,74. 50. +r, 74.60. Ge 
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In cuprate superconductors with large coherence- 
length anisotropies (7 > 10 2 ), the Josephson coupling 
energy between adjacent layers (or bilayers) is quite 
small. The Josephson plasma frequency lu p , which 
measures the coupling between layers, has been ob- 
served in c-axis reflectivity experiments at far-infrared 
frequencies (uj p /2it ~ 30 — 70cm~ 1 in the single-layer 
cuprate La2_ x Sr x Cu04) 0. In the bilayer cuprate 
Bi 2 Sr 2 CaCu208 (Bi 2212), the large separation s be- 
tween adjacent bilayers (s ~ 15. 8 A) weakens the Joseph- 
son coupling significantly. An applied magnetic field H 
further suppresses uj p to frequencies in the microwave 
range (10-100 GHz; 1 cm" 1 ~ 30 GHz). The plasma 
resonance in Bi 2212 is readily observed as sharp reso- 
nances in the absorption of incident microwave power as 
H is slowly swept . Empirically, the field dependence 
of uj p may be described by lu 2 ~ H~^, with // ~ 0.8 when 
H||cg|§. 

Recently, the identification of the resonant absorption 
with excitation of the plasma mode has been called into 
question by Sonin || , who discusses a number of difficul- 
ties with the plasma mode interpretation P~p|,p|,^0t . To 
help distinguish between competing models, quantitative 
comparisons between experiment and theory seem desir- 
able. The strong variation of the resonance field with 
field tilt angle 9 (the angle between H and the ab plane) 
is a highly specific feature of this phenomenon. We re- 
port measurements over a wide range of fields and 9, and 
compare the data with a recent model by Koshelev |ll[ . 

We briefly review the Josephson plasma mode. The 
transfer of n p Cooper pairs between two Josephson- 
coupled superconducting layers leads to a charging en- 
ergy cost of (2n p e) 2 /2C, where e is the elementary 
charge, and C the unit-area capacitance. In addition, 
deviation of the relative phase Sip from zero incurs a cost 
in energy Ej(l — cosSip) where Ej is the Josephson en- 



ergy (per unit area). Hence, the Hamiltonian with n p 
and Sip as conjugate variables is (for H = 0) |l^] 

Hj = (2n p e) 2 /2C + Ej{1 -cos Sip). (1) 

In analogy with the pendulum, the charging energy oscil- 
lates out-of-phase with the Josephson energy at the reso- 
nance frequency ui p given by uj p — (2tt / 4>q) 2 E j / C where 
4>q = h/2e is the flux quantum (h is Planck's constant). 
In finite field, pancake vortices in the layers strongly dis- 
tort the local phase by adding a component ip n (r) that is 
static on the time scale of lu p (where r = (x. y) and n is 
the layer index; we choose axes with c || z and H x c || x). 
If the vortex array is highly disordered, ip n (r) fluctuates 
stochastically with r and n. Adding this static phase to 
the dynamic phase variable in Eq. [j], we find that w 2 is 
reduced p|JlO|| by the factor (cos<^ n)n+ i) (thermally aver- 
aged and integrated over the sample) , where 



0«,n+i(r) = ¥>n+i(r) - p>n(r) - (27t/</> ) 



di-A (2) 



is the gauge-invariant phase difference between layers n 
and n+ 1 with dt || c, and A is the vector potential. The 
factor (cos 4> n ^ n +i) contains the main effect of field on 
Lu p . In general, quantitative comparisons of calculations 
1) with experiment are difficult because they 
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require detailed knowledge of vortex configurations. 

If the Josephson energy is small, however, the vortex 
positions at high temperatures are very weakly corre- 
lated between adjacent layers. In the presence of a tilted 
field, the flux density parallel to the layers is nearly uni- 
form and close to the applied in-plane field induction 
By. While the vortices in each layer produced by B z 
remain highly disordered, the in-plane component cou- 
ples to 0„ in+ i uniformly throughout the sample. This 
situation provides a more transparent comparison with 
experiment. 
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First, we assume H || c. In the high-temperature 
limit, complete randomness in vortex positions should 
suppress (cos 4> n ,n+i) to a value close to zero. As T 
decreases, we expect correlations of the vortex posi- 
tions to increase (within each layer and between layers). 
Clearly, the magnitude of (cos ) grows in propor- 

tion to the strength of these correlations. The relation 
between (cos 4> n ,n+i) an d the phase correlation function 
(exp[i(0„ >n+ i(r) - </V„ + i(0))]) o is derived by expand- 
ing the partition function in the small parameter f3Eja 2 
(where = l/k B T and a 2 = 4>q/B z ) gl|. We find 

(cos0„ !T!+ i) ~ — / d 2 r(exp[i(0 n! „ + i(r) - ^„, n +i(0))]} . 

(3) 

As to 2 ~ (cos</>„ iIl+ i), Eq.|| states that the plasma fre- 
quency (squared) is the average of the phase correlation 
function over the sample. By the definition of <fi n .n+i, 
it is clear that the correlation involves <p n at 4 points (2 
points separated by r in each layer). If in-plane correla- 
tion dominates the interlayer correlation, the integrand 
in Eq.|| reduces to the modulus-squared of the single- 
layer correlation function S\{v) — (exp[z(</?(r) — y>(0))])o 

0- 

Next, we consider tilting the field. If B y is assumed 
to be uniform, it simply adds to the phase difference in 
Eq.|^ a term proportional to the flux threading the area 
defined by the 4 points, viz. {2tt / 4>o)B y sx [the A term 
in Eq. 0] . Equation then simplifies to 



plane, while the large-ga range is accessed with an intense 
H with 8 -v 0. 

Finally, with the approximation |Si(r)| 2 ~ exp(— r 2 /£ 2 ), 
we have Koshelev's result InTl 



(COS (p n ,n+l 



PEj 



d 2 r\Sx{r)\ 



(4) 



where q = 2irsB y /tpQ. We have the key result that B y 
adds a phase shift that increases linearly with x with 
the wavevector q. Upon integrating over the volume, 
we obtain the Fourier transform of |Si(r)| 2 , which may 
be observed as uj 2 versus B y . Thus, in a tilted field, 
information on the in-plane correlation of the vortices, 
inserted by B z , is obtained by measuring u> p versus B y . 

It is instructive to compare this situation with the be- 
havior of an ideal, low-T c junction in a field. There, 
the Josephson coupling is uniform over the length L of 
the junction. The field produces a Fresnel interference 
pattern with period fixed by L, identical to that pro- 
duced by the diffraction of light through a slit of uniform 
transparency. In the present case, however, vortex disor- 
der strongly modulates the Josephson coupling on some 
correlation length-scale. The Fourier pattern produced 
is closer in analogy to Bragg diffraction from a liquid. 
The predicted behavior of uj 2 versus 1 /6 is analogous to 
the variation of the structure factor <S(q) versus q mea- 
sured in a liquid with strong correlation. From the re- 
lation qa = 2ns(B y /(j)o tan^) 1 / 2 , we find that the small 
qa regime is explored in a weak H tilted away from the 
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where £ is the in-plane correlation length, and / = 
a -2 J d 2 r\Si(r)\ 2 — 7r£ 2 /a 2 equals the number of vortex 
pancakes within a correlated patch. We compare Eq.|5| 
with our experiment. 

Crystals of Bi2Sr2CaCu20g with transition tempera- 
ture T c ~ 92 K were cut from a boule grown in an im- 
age furnace and annealed in flowing oxygen to obtain 
optimum oxygen content. The samples were glued to a 
diamond substrate and placed in an oversized circular 
waveguide. Excitation of the plasma mode was detected 
bolometrically as sharp resonances in the absorption of 
microwave power versus field. The absorption was de- 
tected by a cernox sensor. 

At each tilt angle 8, the resonance field Bq is deter- 
mined by slowly sweeping the field at fixed microwave 
frequency uj The inset in Fig. [I] displays sweep-up 

and sweep-down traces at 50 K and 48 GHz, for three 
values of 8. By repeating the sweep at 4 to 6 values of ui 
(between 30 and 70 GHz), we map out the dependence 
of ojp on B at that 8. The procedure is repeated at ~ 20 
values of 8 to generate the full dependence of Bq on u> 
and 8. The main panel of Fig. [I] shows the variation of 
B with 8 at five values of ui. With decreasing 8, Bq 
rises to a maximum value, and then decreases linearly to 
a cusp at alignment. In Fig.|^, we display a different cut 
of the data by plotting u> versus Bq at fixed 8 in log-log 
plot. The solid lines in both figures are fits to a slightly 
modified form of Koshelev's expression which we discuss 
below. 

We summarize the key features observed. Over our 
range of cj, it is convenient to approximate the field vari- 
ation of Up by the power-law, lo 2 = AB^ 11 , where A and /i 
are empirical parameters that depend on 8 and T. When 
H is normal to the layers, the power-law exponent of uj 2 
vs H is fractional, viz. fi ~ 0.8. For small tilts of H away 
from the c-axis, the power law with the same fi and A 
remains valid, if we replace H with H z , as expected for a 
highly anisotropic system. However, as H is tilted closer 
to the plane {8 < 15°), strong deviations from this sim- 
ple scaling become apparent. While the data may still be 
described by a power law, the exponent \i rises rapidly 
to 1.80 as 8 — > 0. The plasma frequency u) p , observed at 
fixed H, increases rapidly to a maximum at 8 near 1°, 
and then decreases linearly with \8\ to a finite value as 
8^0. For experimental reasons, it is convenient to plot 
Bq versus 8 at fixed to. Nonetheless, the same trend is 
observed. Bq rises to a maximum at 6 ~ 1°, and then 
decreases to a finite value at 8 = 0. In this narrow range 
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of 8, the empirical parameters \x and A change rapidly, 
reflecting the increasing influence of the Josephson vor- 
tices. 

As it stands, Koshelev's expression Eq.|| is in close 
agreement with the observed dependence of Bq on 0, if 
to is kept fixed while 8 is varied (because u> is the inde- 
pendent experimental variable, we need to invert Eq. ^]to 
calculate the B versus 8 profile). Fits to the measured 
Bq at angles 8 > 1° are shown as solid lines in the main 
panel of Fig. [I] using Ej and / as adjustable parameters. 
However, if we vary cj, in addition to 8, the original Eq.|5| 
is inadequate. We find that / should have a field depen- 
dence of the form 



f(B)=g(Bsmt 



(0) 



with a a small, positive exponent and g a constant. 

This additional B dependence is immediately apparent 
in the limit 8 = 90° (H || c). Our measurements show 
that uip varies as B~ 8 , whereas Eq. || predicts to 2 ~ 
Ejfa 2 /(2k B T). This implies that / - B° z - 2 , as in Eq.§ 
with a ~ 0.2. 

The finite value of the anomalous exponent a is ap- 
parent at other tilt angles as well. To bring this out, 
it is preferable to plot, in place of oj 2 , the quantity 
W(B,8) = (oj p /2ir) 2 B sin 6*, in which the trivial depen- 
dence on B z is removed. We show in Fig. |] fits of W(B, 8) 
to the expression 



W(B,8) = M(B sin (9)" exp 



^(.Bsin^Trs^cosfl) 2 



Do sin f 



(7) 



The fitted data (54 points) are restricted to angles be- 
tween 0.8° < 8 < 24°, and frequencies between 40 and 
60 GHz. From a least-squares fit, we find at 30 K, a 
= 0.25, M = 1.31 xlO 20 , and g = 1.52 (M and g are 
in SI units). The corresponding numbers at 50 K are 
0.19, 6.52 xlO 20 , and 1.44. We remark that the depen- 
dence of W on B and 8 is non-trivial. Nonetheless, the 
expression Eq. ^ successfully describes its behavior with 
only 3 adjustable constants at each T. (The curves would 
be simple exponential decays if a were zero). Consider- 
ing the large range of 8 covered, we conclude that the 
fit provides strong support for the modified Koshelev ex- 
pression, with f(B) given by Eq.||. 

The quantity f(B) is the number of vortices inside a 
correlated area. The experiment shows that, as B z in- 
creases, f(B) increases slightly faster than a -2 , which 
implies that the number of vortices inside a correlated 
area increases slowly with field. This is in accord with 
the physical expectation that correlation is relatively 
stronger at denser packing. The fits provide a quanti- 
tative description of the vortex liquid state at 30 and 
50 K. In terms of M and g, the Josephson energy is 
given by E 2 = 2fcsT(M0o/ff)( £ o£r/s), where £q and e r 



are the vacuum permittivity and dielectric constant, re- 
spectively. With e r ~ 25, we obtain the value Ej = 1.50 
and 1.37 xl0~ 8 J/m 2 at 30 and 50 K, respectively. These 
values correspond to zero-field critical current densities 
J c0 = 2nEj/(t)v = 4,600 and 4,160 A/cm 2 . Between 50 
and 30 K, the small parameter (3Eja 2 nearly doubles 
from 0.041 to 0.075 (in a perpendicular field B z = 1 T). 
Nevertheless, the perturbation expansion remains valid. 
With / given by gB z , we calculate (cos <t> n ,n+i) = 0.050 
and 0.028 at 30 and 50 K, respectively, in a field of 1 T. 

The picture that emerges is that, at B z = 1 Tesla, an 
average of 1.5 (1.4) pancakes is captured in a correlated 
patch at 30 (50) K. These numbers increase slowly with 
field as B a , implying that the correlation length mea- 
sured in units of a increases slowly as the packing den- 
sity in the liquid increases. The smallncss of (cos <p n ,n+i) 
confirms that adjacent layers are only weakly coupled. In 
a perpendicular field, this number is further reduced as 
a power law B~ - 8 . 

We briefly discuss the large-g (small distance) limit. As 
is apparent in Fig. |l|, the fits deviate significantly from the 
measurements below 1°. In Eq. ^, ui p goes to zero rapidly 
as 8 — > 0, whereas the measured Bq decreases linearly to 
a finite value. Using our empirical power-law description, 
we may express u> 2 as B y 2 ^ , where the exponent r\ ~ 
0.1 is again fractional. Our current understanding is that 
at short distances (r < £), the Gaussian approximation 
for Si is not valid. The correlation function seems to 
increase much faster than in the Gaussian expression. 
This point requires more investigation. 

In summary, we have carried out a detailed investi- 
gation of how the Josephson plasma resonance field Bq 
varies with tilt angle and microwave frequency at 30 and 
50 K. At fixed to, Koshelev's expression Oj describes the 
angular dependence remarkably well for 8 larger than 1°. 
By introducing a weak field dependence to the quantity /, 
which counts the number of pancakes inside a correlated 
patch, we are also able to account for the field dependence 
of Up. In this experiment a large number of independent 
measurements have been used to pin down the parame- 
ters / and Ej. Thus, it yields a rather reliable estimate 
of these quantities, which characterize, respectively, the 
intraplane and interplane behavior of the vortex liquid 
state. Our results also relate quite naturally the anoma- 
lous exponent 1 — \i to the growth of the correlated area 
with B z . 
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FIG. 1. (Main panel) The dependence of the resonance field 
Bo on tilt angle 9 and microwave frequency uj in Bi 2212 
at 50 K (9 is the angle between H and the ab plane). At 
each value of 9 and ui, Bo is determined by slowly scanning 
H. As 9 decreases, Bo rises to a maximum near 1°. Solid 
lines are fits of data above 1° to Eq.|^. In the 'scattering' 
language, the region 9 > (<) 1° corresponds to small (large) 
q£. Broken lines are continuations of fits to angles below 1°. 
The inset shows absorption curves versus H at 50 K at the 
angles indicated (uj/2n = 48 GHz). Arrows indicate scan 
direction. 

FIG. 2. The field dependence of the plasma resonance fre- 
quency in Bi 2212 at selected values of 9 measured at 30 K. 
Solid lines are fits to Eq.|^. Data at 50 K are similar. 

FIG. 3. The field dependence of W = {uj p /2ty) 2 B sin 9 at 
selected values of tilt angle at 30 K (upper panel) and 50 K 
(lower). The angular dependence above 1° in Fig. 1 is domi- 
nated by (Bo sin#) -1 . By removing this factor in W, we iso- 
late the residual effects of the anomalous exponent a (all the 
curves would decrease exponentially with Bo if a = 0). The 
solid lines are fits to Eq.|. At 30 K (50 K), a = 0.25 (0.19), 
M = 1.31 (6.52) x 10 20 , and g = 1.52 (1.44) in SI units. 
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